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TECHNICAL MEMORANDUM 1239 


TWO-DIMENSIONAL MOTION OF A GAS AT LARGE SUPERSONIC VELOCITIES* 

By S. V. Falkovich 


A large number of papers have been devoted to the problem of 
integration of equations of two-dimensional steady nonvortical 
adiabatic motion of a gas. Most of these papers are based on the 
application of the hodograph method of S. A. Chaplygin in which 
the plane of the hodograph of the velocity is taken as the region 
of variation of the independent variables in the equations of 
motion; the equations become linear in this plane. The exact 
integration of these equations is, however, obtained in the form 
of infinite series containing hypergeometric functions. The 
obtaining of such solutions and their investigation involve exten- 
sive computations. As a result, methods have been developed for 
the approximate integration of the equations of motion first 
transformed to a linear form. S. A. Chaplygin in reference 1 
first pointed out such an approximate method applicable to flows 
in which the Mach number does not exceed 0.4. 

S. A. Christianovich (reference 2), in solving the problem 
of the flow with circulation about a wing in a supersonic stream, 
gave as a first approximation a generalization of the method of 
Chaplygin to the case where the region of variation of the velocity 
in the hodograph plane lies within a sufficiently narrow ring 
entirely inside the circle W < a* . At the same time and inde- 
pendently an analogous method was proposed by Tsien and von Karmari 
(references 3 and 4). These methods are not applicable for 
Mach numbers near unity. In the papers by F. I. Frankl (refer- 
ence 5) and S. V. Falkovich (reference 6), approximate equations 
of motion are given suitable for the investigation of flows in 
passing through the velocity of sound. In his recently published 
paper, S. A. Christianovich (reference 7) showed that for Mach 
numbers 1.05 <M <3.5 the equations for the stream function and 
the velocity potential may approximately be replaced by the equa- 
tions of Darboux with integral coefficients and hence the general 
integral obtained in finite form. 


*"Institut Mekhaniki Akademii Nauk Souiza, SSR" Prikladnaia 
Matematika I Mekhanika, Tom XI, 1947. 
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The equation of motion of a gas is investigated herein for 
large supersonic velocities and a method is shown for the approxi- 
mate integration, which gives sufficient accuracy for Mach numbers 
greater than 4. 

1. Fundamental equation and its transformation. - The equa- 
tion determining the velocity potential of the two-dimensional 
irrotational adiabatic motion of a gas has, as known, the form 

(1 ’ 1) 


To this equation the transformation of Legendre is applied. 

As the independent variables in place of x and y are taken 
Bcp/Bx = u and Bcp/By = v and in place of the velocity poten- 
tial qp(x,y), the Legendre potential is introduced. 

$ (u,v) = ux + vy - cp(x,y) (1.2) 

Equation (1.1) becomes 

( “ 2 ‘ t2) 0 + 2UT ^ + (a2 ' “ 2> 0 - 0 (1 - S) 


If the function $(u,v) is determined, then by differen- 
tiation of equation (1.2), the coordinates x and y of the 
corresponding point in the flow plane are determined. 


x = BG/Bu y = B$/Bv (1.4) 


Pass in the hodograph plane u,v to polar coordinates and 
set u = W cos 9 and v = W sin 9 } equations (1.3) and (1.4) 
then become 


a 2 - W 2 B 2 $ B^t a 2 - W 2 B4> = Q 

a 2 W 2 Be 2 BW 2 a 2 W BW 


(1.5) 


x 


cos 



sin 6 B$> 

w Se 


7 


sin 9 



cos 0 5$ 

W Se 


( 1 . 6 ) 


If only supersonic velocities are considered, it is more con- 
venient in place of the independent variable W to introduce in 
equation (1.5) the new variable 



1185 


HACA TM 1239 


3 


z = 1 /a/m 2 - 1 


(1.7) 


Thus z will be small for large Mach numbers and will increase 
indefinitely as M— >1. 

In order to carry out this transformation, from the equation 
of Bernoulli the expression of W is found in terms of z. Thus 


h 2 a* 


( h 2 . |±i = 6 ) 


1 1 h2 % 2 

2 (K-l) M 2 2W 2 


If M 2 is replaced by z, from equation (1.7) 

& W) , 

w 2 ? ^ 

1+lrz 2 

By making use of this expression, pass in equations (1.5) and 
(1.6) from the variable W to the variable z; then 

(h 2 -l) 2 9 2 4> . 3h 2 z 4 + 2h 2 z 2 + h 2 - 2 


(h 2 -!) 2 d 2 $ + 5h 2 z 4 + 2h 2 z 2 + h 2 - 2 9$ = Q 

9z 2 (l+z 2 ) 2 (l+h 2 z 2 ) 2 be 2 (1+z 2 ) (l+h 2 z 2 )z 9z 


x.-i h±£Ui*!£l C08 9 ft + Ein a n 

w (h 2 -i)z 9z 9e 


-it 

“ ' w 


<il9 L». oos e 

(h -l) z oz 


The characteristics of equation (1.9) may be taken in the form 
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x = e 



(h 2 -l)dz 
(l+z 2 ) (l+h 2 z 2 ) 


4 



(h-l)dz 

(l+z 2 ) (l+h 2 z 2 ) 


The line of maximum velocity z = 0 in the hodograph plane 
will correspond in the plane of the characteristics Xu to the 
line u - X = 0. If the integration is carried out, 

X * 0 - (h arc tg hz - arc tg z) p = 9 + (h arc tg hz - arc tg 

( 1 . 11 ) 


2. Investigation of equation (1.9). - If equation (1.9) is 
referred to the characteristic coordinates X u> 


d 2 $ h 2 - 2 - 2z 2 - i^z^ /S $ _ S gA _ n 

5A 4 (h 2 - l)z 


( 2 . 1 ) 


From equations (1.11), 

U - X = 2(h arc tg hz - arc tg z) (2.2) 


From equation (2.2) it follows that the coefficient of the 
equation (2.1) 


-r / h 2 - 2 - 2z 2 - h 2 z 4 „ \ 

— 7 (S ' - ~ f)W (2 - 3) 

is a function of the difference u - X. Equations (2.2) and (2.3) 
give this dependence in parametric form. 

From equation (2.3) it follows that the function L(u - X) 
is negative for z 2 > (h 2 - 2)/h 2 and positive for 

z 2 < (h 2 - 2)/h 2 . For z 2 = (h 2 - 2)/h 2 the function L(u - X) 
becomes zero. According to equation (1.7) this function corre- 
sponds to the Mach number 


M = 2/ = 1.565 

S. A. Christianovich showed (reference 8) that for_a given 
value of the Mach number there is a change in the direction of 
curvature of the characteristics in the flow plane. The graph of 
the function L(p - X) is shown in figure 1. 
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If the right side of equation (2.2) is expanded in a series 
in powers of z, 

E = (k 2 - 1) z - z 3 + . . . (2.4) 


from which 


z 


u - A + h 2 + 1 
2(h 2 - 1) 24(h 2 - l ) 3 


(r - A ) 3 + 


If this series is substituted in equation (2.3), 


L(u - A) 


h 2 - 2 _ (h 2 + 1 ) (h 2 - 2 ) + 6 

2 (^ A) 24 (h 2 - l ) 2 


(|x - A) + . . . 


Set h 2 = 6 . Then 


L(jjl - A) = — ~ - 0.057 (p - A) + . . . 

{-I — A 


For M > 4, assume 


L(ti - A) * 5 -^ 


Equation (2.1) then becomes 

d 2 $ _ 2 SS>\ _ _ 

aX Sp X - p ^p/ 


(2.5) 


Equation (2.5) is the equation of Darboux with integral 
coefficient. The general integral of this equation has the form 


5 2 X(A) - Y(u) = X«(A) + Y'(p) _ X(A) - Y(p) 
' SA Sp A - V- (A - p ) 2 (A - p ) 3 


( 2 . 6 ) 

where X(A) and Y(p) are arbitrary functions of their arguments. 
The expression (2.6) is the asymptotic integral of the exact 
equation ( 2 . 1 ) for z— > 0 , that is, for M -*•<». 
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Equations (1.10) for the coordinates x and y. after 

becomes on the 

-A 


w 


y = + i 


'I® + 8 i n h + u . 1 /a® a$\ 

2 2 ( " £a/ cos 


a_+_ 

2 


^ + |®Y OB ?L±Ji _ 1 fa d$\. X + u 

2 ^U'W sin z 


> (2.7) 


From equation (2.2) for small values of z. 


J 


z = — ^ ^ = 4 ~ A 

2(h 2 - 1) 10 

If v?S ±S value of 2 is substituted in equation (2 7 ) and 

SS^S y fT’i 1 °“ from eqieMon (2 ?) 

»ith the aid of equation (2.6) the final solutii is “btai^i. 


m = .i 
w 


^r_+_r _ x« - yA a + u 
x(x-n) 2 (x-u) 3 / 1 ~T M 




7 ' + i 


XI + Y" _ 2 x» - y» N 
V(^) 2 (A-n) 3 / 


cos + 


10 


X - Y" _ x* + T . y .. Y\ * 

\(*-n) 3 " ( A -n )4 + 12 7^)5) 8in — 


and ^rSm^ncX: 1 ^^ 9 VelOClty »*-«- <PCw,0) 

^ = £2 (m2 ,1 1W 

dd p Sw aw p ^ w Se 


♦ 


lies 
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If the variable z is introduced in place of the modulus of the 
velocity W according to equation (1.7) and equation (1.8) used, 
these equations are transformed to the form 

K K k+1 "\ 

^ - (1 . h 2 * 2 )*- 1 (S 2 - I) 1 -" (1 + z 2 ) F* H 


2-K K+1 

|| - - (1 + hV)'* 1 (1 + Z 2 )' 1 (h 2 - l)" 1 z 1 ^ || 


>( 2 . 8 ) 




For small z, equations (2.8) may be replaced approximately by 
the equations 

K K+1 

^2 = - (h 2 - l ) 1- K z 1-K ^ 

58 2 5z 


= - (h 2 - 1) 


2-K Kfl 
K-l 1-K 5ty 


Si 


If the velocity potential cp is eliminated, the equation for the 
stream function \J/ is 


(t? - D 2 - 2f| + £ |s . o 

■ " ° z 5z 


56* 5z 


(2.9) 


the characteristics of which have the form 


X - 0 - (h 2 - l)z \i m 0 + (h 2 - l)z 


If equation (2.9) is referred to the characteristics, it is trans- 
formed into the form 


+ it > 2 , . 0 
5X 5jjl 2(X - u) \oX ayi J 

Setting h2 a 6 gives 
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The equation of Darboux with integral coefficient, which is anal- 
ogous to equation (2.5), is integrated in finite form. 


5. Criterions of similarity. - By examining equation (1.9) 
for Mach numbers near 1 (z— and also for large Mach n umb ers 

(z— »0), certain criterions of similarity can be established 
that may be useful in evaluating experimental data obtained in 
wind-tunnel tests. For large values of z, equation (1.9) can 
be replaced by the approximate equation 


a 2 $ _ (h 2 - D 2 a 2 $ , 5 50 _ Q 
Sz 2 ’ h 4 z 8 de 2 * = 


(3.1) 


and for small values of z by the approximate equation 


d 2 $ 

dz 2 



+ h 2 - 2 

be 2 z 'Sz 


0 


(3.2) 


A thin slightly cambered airfoil at small angle of attack is 
now considered in a plane-parallel nonvertical gas flow with Mach 
number Mq at a large distance from the airfoil. The profile 

chord is denoted by Z and its maximum thickness is denoted by 6. 
If it is assumed that is near unity, in equation (3.1) 

9 = 6*b/l z = z*z Q ^z 0 = 1/y^Tl) 

and equation (3.1) becomes 


where 


rr 2 h 2 - 1 5 M , 

1 ST 5 h^z* 8 ae* 2 ***^ 


% 3 | /(M 2 - l) 3 


(3.3) 


(3.4) 


The number K ^ may be called the criterion of similarity in 

the sense that two flows having different Mach numbers about two 
airfoils of different thicknesses 6 and different chords Z but 
for which the values of are the same will be determined by 

the same '* nondimensional 11 equation (3.3); hence, for the drag 
coefficient 


- t (%) 
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In an entirely analogous manner starting from equation (3.2), 
a second criterion of similarity valid for large Mach numbers is 


*2 - f A> 2 - 1 s H 

This criterion of similarity has recently been obtained by Tsien 
by a different considerably more complicated method. 


Translated by Samuel Reiss 

National Advisory Committee 

for Aeronautics i 
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